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Abstract 

We use dynamic light scattering and computer simulations to study equilibrium dynamics and 
dynamic heterogeneity in concentrated suspensions of colloidal hard spheres. Our study covers an 
unprecedented density range and spans seven decades in structural relaxation time, r a , including 
equilibrium measurements above ip c , the location of the glass transition deduced from fitting our 
data to mode-coupling theory. Instead of falling out of equilibrium, the system remains ergodic 
above (p c and enters a new dynamical regime where r a increases with a functional form that was not 
anticipated by previous experiments, while the amplitude of dynamic heterogeneity grows slower 
than a power law with r a , as found in molecular glass- formers close to the glass transition. 

PACS numbers: 05.10.-a, 05.20.Jj, 64.70.P- 
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Hard sphere assemblies often constitute the simplest model to tackle a variety of fun- 
damental questions in science, from phase transitions in condensed matter physics to the 
mathematics of packing or optimization problems in computer science. Experimentally, hard 
spheres systems are obtained using colloidal particles [l| , emulsions, or granular materials . 
When crystallization is avoided, e.g. by introducing size polydispersity, hard spheres at ther- 
mal equilibrium become very viscous and eventually form an amorphous solid 3| at large 
volume fraction, ip, in analogy to the glass transition of molecular liquids {4] and the jam- 
ming transition of grains {2]. However, the nature of the colloidal glass transition, its precise 
location, the functional form of the structural relaxation time divergence, and the connection 
between slow dynamics and kinetic heterogeneities remain largely open issues 5|, |6j. 

For hard spheres at thermal equilibrium, several distinct glass transition scenarios have 
been described. In the first, the viscosity or, equivalently, the timescale for structural 
relaxation, T a (cp), diverges algebraically: 



(1) 



This is predicted [7J by mode coupling theory (MCT), and supported by the largest set of 
light scattering data to date fjj]. Packing fractions (p c ~ 0.57 — 0.59 are the most often quoted 
values for the location of the 'colloidal glass transition'. It is widely believed that a truly 
non-ergodic state is obtained at larger (p . Within MCT, the amplitude of dynamic 

heterogeneity quantified by multi-point correlation functions also diverges algebraically. In 
particular, the four-point dynamic susceptibility should diverge as Q|: Xa ~ (v^c ~ i P)~ 2 ~ 
t^ 7 , a prediction that has not been tested experimentally. 



Several alternative scenarios 

T a {<p) = ^exp 



11] suggest a stronger divergence: 
A 



(2) 



Equation (jSJ) with 5 = 1 is frequently used to account for viscosity data {(| because it 
resembles the Vogel-Fulcher-Tammann (VFT) form used to fit the viscosity of molecular 
glass- formers {4], with temperature replaced by <p. Moreover, it is theoretically expected 
on the basis of free volume arguments 9], which lead to the identification ipo = (f rcp , the 
random close packing fraction where osmotic pressure diverges. Kinetic arrest must occur 
at ip rrv (possibly with 5 = 2 [u]]), because all particles block each other at that density 10, 



12 



2 [10j), because all particles b 
131 ]. Entropy-based theories and replica calculations 



111 ] predict instead a divergence of 



r a at an ideal glass transition at (fo < (f vcp , where the configurational entropy vanishes but 
the pressure is still finite. Here, the connection to dynamical properties is made through 
nucleation arguments [ljj] yielding Eq. (j2J), with 6 not necessarily equal to unity In this 
context, the amplitude of dynamic heterogeneity should increase only moderately, typically 
logarithmically slowly in r a fl6 |. 

In molecular glass-formers where dynamical slowing down can be followed over as many 
as 15 decades, the transition from an MCT regime, Eq. (JTJ, to an activated one, Eq. (j2J), 
has been experimentally demonstrated \&\. For colloidal hard spheres, the situation remains 
controversial, because dynamic data are available over a much smaller range 
typically five decades or less. Crucially, equilibrium measurements were reported only for 
<p < (p c , leaving unknown the precise nature and location of the divergence. Theoretical 
claims exist that the cutoff mechanism suppressing the MCT divergence in molecular systems 
is inefficient in colloids due to the Brownian nature of the microscopic dynamics, suggesting 
that MCT could be virtually exact 18[. This viewpoint is challenged by more recent MCT 
calculations 19j, and by computer studies of simple model systems where MCT transitions 
are avoided both for stochastic and Newtonian dynamics 20l. |21|. 

Here, we settle several of the above issues by studying the equilibrium dynamics of col- 
loidal hard spheres using dynamic light scattering and computer simulations. By extending 
previous data by at least two orders of magnitude in r a , we establish that the volume frac- 
tion dependence of both r a and %4 follows MCT predictions only in a restricted density 
range below our fitted ip c ~ 0.59. Unlike previous studies, we provide equilibrium measure- 
ments above tp c , thereby proving unambiguously that in our sample the algebraic divergence 
at (f c is absent. Instead, a new regime is entered at larger </?, where the dynamics is well 
described by Eq. (j2J) with 5 ~ 2 and <po much larger than tp c . The amplitude of kinetic 
heterogeneities then grows slower than a power law with r a , as in molecular glasses close to 
the glass transition. 

Dynamic light scattering (DLS) experiments are performed in the range 0.01 < cp < 
0.5981. We use poly-(methyl methacrylate) (PMMA) particles of average diameter o = 
260 nm, stabilized by a thin layer of grafted poly-(12-hydroxy stearic acid) (PHSA). The 
size polydispersity, about 10%, is large enough to prevent crystallization on a timescale of 
at least several months. The particles are suspended in a mixture of cis-decalin and tetralin 
that almost perfectly matches their average refractive index, allowing the dynamics to be 
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probed by DLS. Additionally, a careful analysis of the combined effects of optical and size 
polydispersity shows that we probe essentially the self-part of the intermediate scattering 
function 221 ] : 

F a {q,t) = (^p^<'^-'^y (3) 

Here, ij(t) is the position of particle j at time t, q is the scattering vector (qa = 6.5, 
close to the first peak of the static structure factor) and brackets indicate an ensemble 



average. A combination of traditional 



23] and multispeckle [24j DLS is used to measure 



the full decay of F s (q,t). We carefully check equilibration by following the evolution of 
the dynamics after initialization, until F s stops changing over a time window of at least 
10r Q . Samples are prepared by dilution, starting from a very concentrated batch obtained 
by centrifugation. All volume fractions relative to that of the initial batch are obtained with 
a relative accuracy of 10~ 4 , using an analytical balance and literature values for particle and 
solvent densities 17] • Relative volume fractions are converted to absolute ones by comparing 



the experimental <p dependence of the s 
to two sets of theoretical calculations 



rort-time self- diffusion coefficient measured by DLS 



251 ] at low density, cp < 0.2. For less polydisperse 



samples, this calibration method yields ip values compatible with those obtained by mapping 
the experimental freezing fraction to ipt = 0.494 26j. The remaining uncertainty on the 
absolute cp is about 5%, because j^J contains two slightly different predictions. To ease the 
comparison with the simulations, we set the absolute tp, within this uncertainty range, so 
that our experimental and numerical r a closely overlap for tp > 0.55. 



We use a standard Monte Carlo algorithm 



211 ] to study numerically a 50:50 binary mixture 



of hard spheres of diameter a and lAa, known to efficiently prevent crystallization. We work 
in a three dimensional space with periodic boundary conditions, and mainly use N = 10 3 
particles. No noticeable finite size effects were found in runs with N = 8 • 10 3 particles 
performed for selected state points. In an elementary move, a particle is chosen at random 
and assigned a random displacement drawn within a cubic box of linear size O.lcx centered 
around the origin. The move is accepted if the hard sphere constraint remains satisfied. 
One Monte Carlo step corresponds to N such attempts. The dynamics is characterized by 
the self- intermediate scattering function, Eq. (J3j), measured for qa = 6.1, close to the first 
diffraction peak. 

Representative F s (q, t) obtained by DLS are plotted in Fig. [H showing that the relaxation 
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FIG. 1: (Color online) Time dependence of the self-intermediate scattering function F s (q, t) in DLS 
experiments at qa = 6.5 for representative volume fractions. Lines are stretched exponential fits to 
the final decay, yielding relaxation times spanning about 7 decades. Ergodicity is preserved above 
the (avoided) MCT glass transition at ip c ~ 0.59. 



is fast and monoexponential at low ip, while a two-step decay is observed when increasing p, 
reflecting the increasingly caged motion of particles in dense suspensions We measure 
the structural relaxation time by fitting the final decay of F s to a stretched exponential, 
F a {q,t) = Bexp[-(t/T a )P]. 

Figure [2^ shows r a (ip) for both experiments and simulations. Time units are adjusted 
to maximize the overlap (~ 5.5 decades) of both data sets at high ip. Our experimental 
data are well fitted by Eq. (JTJ in the range 0.49 < tp < 0.585, with tp c = 0.590 ± 0.005 
and 7 = 2.5 ± 0.1. For the slightly less polydisperse sample of Ref. js], a similar power law 
behavior with 7 m 2.7 and <p c = 0.571 — 0.595 was reported, the two quoted values of <p c 
stemming from experimental uncertainty in the volume fraction determination. However, 
our measurements for the largest densities strongly deviate from the MCT fit. Attempts 
to include points at tp > 0.59 in the MCT fit yield unphysically large values of 7. Similar 
deviations are found in our simulations, showing that hydrodynamic interactions play little 
role in experiments performed at large ip, although they probably explain the discrepancy 
with simulations at low volume fraction, see Fig. EK Therefore, our results unambiguously 
demonstrate that the mode-coupling singularity is absent in our hard sphere colloidal system, 



as is also found in molecular glass-formers 



4]- 



What is the fate of the fluid phase above ip c l Figures [2^l and [2fc show that the increase 
of r a at high ip is extremely well described by an exponential divergence, Eq. (J2J). We 
find that the data can be fitted well using the conventional form with 5 — 1, yielding 
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FIG. 2: (Color online) a) Relaxation timescale r a for hard spheres in experiments (black circles) 
and simulations (open triangles), respectively in units of tq = 1 sec and r = 7 - 10 4 MC steps. The 
red dashed line is a power law fit, Eq. dl|), with ip c = 0.590 (vertical dotted line) and 7 = 2.5 ± 0.1. 
The continuous blue line is a fit to DLS data using Eq. ([2]), with cpo = 0.637 and 5 = 2. The zoom 
in the inset shows that the MCT singularity is absent, b) Same data plotted against 1/(1 — ip/ip c ). 
A straight line with slope 7 is obtained in an MCT regime covering almost 3 decades in r a . c) 
Data for (p > 0.41 plotted using reduced variables with ipo = 0.637 and 0.641 for experiments and 
simulations, respectively. 

<Po($ — 1) ~ 0.614 ± 0.002. This is consistent with previous analysis of viscosity data Q|. 
However, the quality of the fit improves when the exponent 5 is allowed to depart from unity. 
The optimal value, robust for both experimental and numerical data, is 5 = 2.0 ±0.2, which 
yields our best estimate for the location of the dynamic glass transition: tpo m 0.637 ± 0.002 
(experiments) and <po ~ 0.641 ±0.002 (simulations). Figure Et shows the linear dependence 
of logr a on ((p — (p)~ 2 , demonstrating the exponential nature of the dynamic singularity. 

The behavior of dynamical heterogeneity provides additional evidence of a crossover from 
a restricted MCT regime to an 'activated' type of dynamics. Using methods detailed in [27, 



2n 



we study the evolution of the three-point dynamic susceptibility defined by: x^(l^) = 
dF s (q,t)/d<p. This linear response function is directly connected to a four-point dynamic 
susceptibility: x^Qit) = N(5F s (q,t) 2 ), where SF s (q,t) denotes the fluctuating part of the 
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self-intermediate function; is a powerful tool to quantify dynamic heterogeneity in glass- 



formers 



161 ]. because it represents the average number of molecules whose dynamics are 



correlated. In hard spheres, the following relation holds 271 ]: 



X4(g, t) = Xa(j9, t)\ v + pk B TK T ((px<p(q, t)) 2 , (4) 

where p is the number density, kt the isothermal compressibility (measured in simulations, 
taken from the Carnahan-Starling equation of state in experiments), and Xi{o.it)\(p denotes 
the value taken by Xi{Q^) m a system where density is strictly fixed. Only the second con- 
tribution to XiiQi t) in (J4]) can be accessed experimentally, but both terms can be determined 



in simulations. We obtain x<p{q^) by applying the chain rule to the fitted F s (B,T a , (3) [28], 
where the p dependence of B,r a ,/3 is fitted by smooth polynomials. Our results are in- 
dependent of the choice of fitting functions, and consistent with that obtained from finite 
differences between data at nearby <p, when available. Figure [3] shows the peak of dynamical 
susceptibilities as a function of tp. First, we numerically establish in Fig. [3h that the term 
comprising x<p is the main contribution to Xi when p > 0.52, implying that a good estimate 
of X4 can be obtained using three-point functions in W spheres, a8 — in 0, and 
established for molecular glass-formers in For both simulations and experiments, the 
MCT prediction for the algebraic divergence of X4(<7, i) only holds over a limited density 
range. Indeed, when plotted against r a , %4 eventually grows slower than a power law, as 
found for the size of dynamically correlated regions in molecular glasses close to the glass 
-a hallmark of activated dynamics jlfij ]. 



transition 



28 



Our results establish the existence of a non-trivial, exponential divergence of T a (tp) at a 
critical volume fraction po m 0.637 much above the putative 'colloidal glass transition' at 
<p c ~ 0.59. It is natural to ask whether tp and p Tcp coincide. This is a difficult question 



because p rcp can always be shifted to a larger value by trading order and packing 12j. For 



the binary mixture studied here, the onset of jamming has been located at p>j = 0.648 13]. 



Furthermore, for 10 % polydispersity, the estimate p Tcp ~ 0.67 was obtained in numerical 



work (29[, well above Pq. Finally, we have employed Monte Carlo simulations to produce 
disordered hard sphere configurations with finite pressure above po by a very fast compres- 
sion of fluid configurations used to produce the equilibrium data in Fig. [1] (open symbols). 



These results support the possibility that p < p rcp , implying a fundamental difference [30] 



between the glass 



11] and jamming [l3j transitions in hard spheres. 
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FIG. 3: (Color online) Peak of dynamic susceptibilites, Eq. (|4|), measured in a) simulations and 
b) experiments. In a) both contributions to \i are compared, validating valuable tool 

to quantify dynamic heterogeneity in hard spheres. The predicted MCT algrebraic divergence 
(red dashed line) holds over a small density range. The inset shows that the size of dynamic 
heterogeneities grows slower than a power law at large r Q , as found in molecular glass-formers. 

In conclusion, we report a set of dynamic data for a well-known colloidal hard sphere sys- 
tem covering an unprecedented dynamic range of equilibrium relaxation timescale. While 
the onset of dynamical slowing can be described by an MCT divergence at a critical volume 
fraction (p c , upon further compression a crossover from an algebraic to an exponential diver- 
gence at a much larger volume fraction <p Q is observed, accompanied by a similar crossover 
for the growth of dynamical correlations. Our results show that the apparent singularity at 
ip c does not correspond to a genuine 'colloidal glass transition', suggesting that the MCT 
transition is generally avoided in colloidal materials, just as in molecular glass- formers. 
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